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Abstract 

In the theory of disordered systems the spectral form factor S(t), 
the Fourier transform of the two-level correlation function with re- 
spect to the difference of energies, is linear for r < r c and constant 
for t > t c . Near zero and near t c its exhibits oscillations which have 
been discussed in several recent papers. In the problems of mesoscopic 
fluctuations and quantum chaos a comparison is often made with ran- 
dom matrix theory. It turns out that, even in the simplest Gaussian 
unitary ensemble, these oscilllations have not yet been studied there. 
For random matrices, the two-level correlation function p(\i,\2) ex- 
hibits several well-known universal properties in the large N limit. 
Its Fourier transform is linear as a consequence of the short distance 
universality of p(\i,\2)- However the cross-over near zero and r c re- 
quires to study these correlations for finite N. For this purpose we use 
an exact contour-integral representation of the two-level correlation 
function which allows us to characterize these cross-over oscillatory 
properties. The method is also extended to the time-dependent case. 
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1 Introduction 



The universal properties of the correlation functions in random matrix the- 
ory have been discussed abundantly and applied to various fields such as 
fluctuations in mesoscopic systems or quantum chaos. 

The Fourier transform of a two-level correlation function is a spectral form 
factor, which we denote by S(t). For the Gaussian unitary ensemble, in the 
large N limit, this form factor S(t) has a simple linear behavior with r up to 
a critical value r c = 2N, beyond which it becomes one [[I]]. This remarkable 
behavior is due to the short distance universality of the two-level correlation 
function p(Ai, A 2 )[@. In the problem of quantum chaos, it is known that the 
level statistics of chaotic systems in a certain energy range, agrees with the 
result of random matrix theory, and the linear behavior of S(t) has been 
derived by the method of perturbation of periodic orbits 0. 

In this article, we evaluate this spectral form factor within random matrix 
theory, in order to characterize the crossover to the linear behavior in the 
large N limit. We will investigate the subdominant term, which has oscilla- 
tory behavior, although it can be neglected in the large N limit. When we 
consider the derivative of S(t) with respect to r, this oscillatory behavior is 
clearly seen, since it becomes of the same order as the linear term. We be- 
lieve that these oscillatory terms near r = and r = 2N, although small, are 
relevant for the discussions of current interest on oscillations in disordered 
metals or in quantum chaos, in the non-universal regions ||, || 

For the discussion of the crossover to the universal linear behavior, we de- 
rive an exact expression for finite N of S(t) . Our analysis is based upon the 
recent calculation of the two- level correlation function [[?[] , in which Kazakov 
contour-integral representation |J has been used. This representation has 
also been used recently for the Laguerre ensemble; it made it easy to char- 
acterize the crossover behavior near the edge and near zero energy for the 
density of state ||. We consider here a similar crossover behavior for the 
two-level correlation function or the spectral form factor. 

We find also, after averaging this form factor over the energy, that the 
corresponding form factor < S(t) > is remarkably related, through a simple 
integration, to the density of state p(r) in the Laguerre ensemble; this density 
is known to posses a universal oscillatory behavior near the origin 10, 11, 

We extend the form factor calculations to the time- dependent case, which 
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is shown to be equivalent to the two-matrix model. In this case, the singular 
behavior at the Heisenberg time is smeared out. 

We further discuss matrix models with an external source for the correla- 
tion functions. We find new characteristic properties of the kernel K N (X,p), 
and its universal behavior. The kernel K^{X,p) has lines of zeros in the 
real (A, p) plane. We briefly study the zeros of the kernels in the two-matrix 
model and in the model with the external source. 



2 Universal behavior of the form factor 

The two-level correlation function p^(X,p) for the random matrix model is 
defined by 

p (2 \X,p) =< ±Tr5(X - M)±Tr8(p - M) > (2.1) 

where the M is an N x N random Hermitian matrix, and the bracket means 
an averaging with respect to the Gaussian distribution; 

P(M) = i eX p(-^TrM 2 ) (2.2) 

The connected correlation function p^(X,p) is obtained by subtraction 
of the disconnected part, which is a product of the density of states p(A) and 
p(/i). This function has a complicated expression with strong oscillations, 
which simplifies only in the short distance limit in which there are a finite 
number of levels between A and p, i. e. for N(X — p) finite in the large N 
limit. Introducing the scaling variable, 

x = 7 ,N{\-p)p{ l -\+ l -p) (2.3) 

and taking the large N limit with a finite x, one finds Jl| 

pf (A, p) ~ h(X - p)p(X) - p(X) P (pf-^ (2.4) 

JX x A 

The spectral form factor S(t) is defined by 

/+oo 
duje iuJT pW(E,E + uj) (2.5) 
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Using the large N, small u limit, we have, leaving aside the delta-function 
term in Q2.4|) , 

^E-^E + ^-^^m ( 2 . 6 ) 



Then, the Fourier integral is evaluated easily, since 

10 , \t\ > 2a 



- e «^) jL;= /i(2a-|t|) ,|t|<2« ^ 



^ Tj "1 ,|r| >2nNp(E) 



This leads to 



Adding the 5-function term of (|2.4j ), we find that S(r) vanishes for r = 0. 
^From this result, we find that if r is order of N, then the integration over 
uj is dominated by a range of order 1/N, and therefore, the approximation 
of p^(X, n) by its short distance behavior (|2.6|) is justified. However, if r is 
order of one, then we have to deal with an integration over a range in which uj 
is not small, and we cannot use anymore the short distance universal behavior 
for p^ 2 ^(A,/i). Therefore, one expects a universal linear behavior in the range 
in which r is of order N. 

We have derived in a previous paper []7| the oscillating short distance 
behavior of ( |2.6| ) by using a method introduced by Kazakov. This method 
gives exact expressions of the correlation functions for finite N. It is very 
convenient for characterizing the crossovers in comparison with the standard 
approach based on orthogonal polynomials. It consists of adding to the 
probability distribution a matrix source, and this external source is set to 
zero at the end of the calculation. (In some cases one is interested in keeping 
a finite external source, as studied recently in 0). We thus modify the 
probability distribution of the matrix by a source A, an iV x iV Hermitian 
matrix with eigenvalues (a 1; ■ ■ • , a^r): 

P A (M) = -^exp(-— TrM 2 - NTiAM) (2.9) 
Z A 2 

We consider the average evolution operator with this modified distribution 

U A (t) =< ^Tre UM > (2.10) 
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The density of state p(X) is its Fourier transform: 

/ + OO rjf 
-oo h e ~ UXUA{t) (2 - n) 

We first integrate over the unitary matrix V which diagonalizes M, (we may 
assume, without loss of generality, that A is a diagonal matrix). This is done 
by the well-known Itzykson-Zuber integral [14||, 



where A (A) is the Van der Monde determinant constructed with the eigen- 
values of A: 

N 

A{A) = l[{a i -a j ) (2.13) 

i<j 

We are then led to 

1 1 N r 

x eM—T.'-' ~ N J2"fi) (2.14) 
After integrating over the r^, we obtain 

U A (t) = 1 £ n f a ~ a ?~h e-&-**- (2.15) 

Instead of summing over N terms, one can write a contour-integral in the 
complex plane, 

1 r rln N v — n — — t 2 
U A (t) = — ff-.m - 1 N )e- Uu ^ (2.16) 

We may now, and only at this stage, let the a 7 go to zero; we obtain 

Uo(t) = --e'^ I %'^{\ - (2.17) 
it J 2m Nu 
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Similarly the two-level correlation function, p^ 2 '(X, //) is obtained from 
the Fourier transform UA{ti,t 2 ), after letting A go to zero 0, 



P (2) (A,/^) = /" f^e"*^^ (2-18) 



where U A (ti, t 2 ) is 



U A {t!,t 2 ) =< I T re^ M ^Tre^ M > (2.19) 



The same procedure leads to 

ufihM) = pE/ n^^ e ^ E( ^ ?+riai)+t(4irai+t2rQ2) (2-20) 



ai,a2 



By integration over the 7*j, we obtain, after subtraction of the disconnected 
part, a representation in terms of an integral over two complex variables 

tt u 4. \ 1 [ dudv _A*L_ iUu _ iUvH ih. N . it 2 , N 

Uoih.h) =--^_e»». -(1 - _)"(!-_) 

1 

X ( M -^-f)( M -,+ f) (2 ' 21) 

where the contours are taken around u = and v — 0. If we let the contour 
include the pole, v = u — %, it gives precisely the disconnnected term U (ti + 
t 2 ), whose Fourier transform is the 5 function part of ( |2 . 4| ) . 

We now write the two-level correlation function as the Fourier transform 
of Uq(\i, X 2 ). In order to show that it takes a factorized from, we shift the 
variables t% and t 2 to t± — > t\ — iuN and t 2 — > t 2 — ivN. Then, one finds 

dh r du it 2 N 1 e _N u 2_g_ it2X2 _ NuXl 
2tt J 2m Nu u + ^ 

= -^K N (X!, X 2 )K N (X 2 , Ax) (2.22) 
We have obtained the integral representation for the kernel K^(X,/j,), 

K N (X, /i) = - f - / —(-^f^re-T (2.23) 
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It may be interesting to note that the same integral expression is obtained 
through the orthogonal polynomial method. In Appendix A, we give this 
derivation. 

The expression of K N (Xi, A 2 ) may be simplified further by the shift t\ — > 

t + ivN, 

K N {\ 1} A 2 ) = / £ f ^1(1 - ttfe-Sr-ivt-^+^-x,) {2M) 

J ZTT J ITTllt Nv 

We may now find the short distance behavior of p^(Ai, A2) in the large 
N limit with a finite value of the variable y = N(X\ — A2). There are several 
procedures to obtain the oscillating universal form. One possibility has been 
discussed in 0. We follow here another procedure for the purpose of later 
use. If we substitute to v, v — > itv, we may then perform the t-integration, 

#*(Ai,A 2 ) = ^/fVr^ ^ 'V-W-) (2-25) 
27r J 2m y ^ — v Nv 

In the large N limit, we may neglect 1/N terms and exponentiate the term 
which is a power of N. We obtain 

We change the contour in the complex plane, and we use the following result, 

/"CO v 

Ki(z) = / e- zcoshx cosh-cix 
2 J 2 

Where Ki(z) is a modified Bessel function. Then, we obtain 



K N (X 1} A 2 ) ~ —e^smiyJ 1 - ^) (2.28) 

The other term Ejv(A 2 , Ai) is obtained in similar way. Thus we get, in the 
large N, finite y, limit, 
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2 2 Sm (^T 



p(Ai,A 2 ) = sin(-5-^ y)sin(-*— ^ y) 



. 2 sin2 (^y^) ( 2 - 29 ) 

We also derive more precise expression for the kernel Kn(X\, A2) from 
^2|). We have 

J 2n J 2m Nv f + 77 

= Nf^i^^e-"f (2.30) 

J 2n J 2ni v + it K ' 



where / is 



t; 2 t 2 



/ = — + — + it\i + v\ 2 - lnt + \nv (2.31) 
The saddle point equations for / become 

df , 1 

(7f V 

df 1 

■77- = t + iX 1 -- = (2.32) 

The four solutions are obtained: v = ie lip , —ie~ tip , t = e~ l9 , —e ld . We define 
Ai = 2sin# and A2 = 2sin<p. The Gaussian fluctuation around the saddle 
point is evaluated by 

1 



v/Scos^ 



e ±qe 



V2cos# 

(2.33) 

Adding these four saddle point contributions, we have 



K N (X 1 ,X 2 ) 



e f(co s2 e-co S 2 V ) cos [ iV (# + ^ + ^ + 2»|£)] cos I(# + ( ^ 



8irNy/cos6cos(p 1 + cos(6 l + <p) 
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+ sin[iV(fl - y + |sin2g - f sin2y)]sinl(g - y) 
1 — cos(# — (f) 

When Ai — A2 is order of 4, we make approximations in fl2.34j ), 

. sin26> sin2o? , 

— <p H — ~ sm(& — (p) + [smu — sin<p)cos0 

~ 7r(Aa - A 2 )p(Aa) (2.35) 

/4— A 2 

where p(Ai) = 2?r 1 , and the denominator of ( |2.34| ) is approximated as 



1 — cos(# — if) ~ |sin 2 (6 l — <p). Then, in the large N limit for fixed iV(Ai — A2) 
we get the short range universal form of ( |2.29| ). Later we will discuss the 
generalization of (|2.34 ) to the time dependent case. 



We now consider the form factor S(r), which is defined by [2.5| . ^From 
the expressions of K N (0,u) and K N (u,0) in ( [2.25| ), we have 

1 f , f dudv e iujT 



8(T) = -J*,f 



In the large N limit, if r is order of N, we may use the previous expressions 
for re-deriving the universal short distance behavior in (2.29p, and obtain 



the linear behavior up to r = 2N . However, for finite N, this function is 
complicated, and we need the study of the oscillating part based on ( |2.36 ). 



3 Oscillatory behavior of the form factor 



We first integrate out u in ( |2.36|) , and by shifting u — > j^u and v — > -kv, we 
obtain 

S(r) = 



dudv _ t 2 

-g 2JVD 



2nN 2 J (2ttz) 2 ^{ Z 2u ) + ( u - 1) 2 (1 - 2v) 



x (l-V(l-V 

U V 



where D is given by 



D 



v 2 (u - l) 2 
1 - 2v + 1 - 2u 



(3.1) 
(3.2) 



9 



A quasi-linear behavior with small oscillations follows from this expres- 
sion. It is interesting first to compute this contour integral ( |3.1|) for finite 
N. In Fig. 1 a) the result is shown for N = 7 . The correction to the linear 
behavior is small, but the derivative of S(t) with respect to r, shown in Fig. 
1 b), becomes of the same order as the constant part, at least near r = 0. 
Returning to the analytic calculation we can obtain exact expressions for this 
oscillatory behavior by a saddle point analysis of (|3.1|) . For this purpose, we 
scale t by r = Nr. Then we have 

S(r) = I f Ud l , 1 e-»' (3.3) 

v^fiV 2 J (2m) 2 ^/ v 2 (l-2u) + (u-l) 2 (l-2v) 

where the exponent / is 

/ = ^ - m(l - -) - ln(l - -) (3.4) 
2D u v 

In the large N limit, we look for the saddle points of u and v in the complex 
plane. They are obtained by 

to = °' to = ° (3 ' 5) 

We get thus the two equations, 

(1 -2uf f_ (\-2vf f_ 

u 2 (l-u) 2 ~ D 2 ' v 2 {l-v) 2 ~ D 2 

As solutions of these equations, we have 

1 - 2u l-2v 



(3.6) 



u(l — u) v(l — v) 



(3.7) 



There are four solutions to these equations, but two of them only are saddle- 
points : a) for the(+) case, u = v. b) for the(-) case, u = T^rf- Although 
v = for the (+) case and v = 1 — u for the (-) case are solutions, they 
are not saddle-points since D vanishes. The case a) u = v leaves us still with 
four different solutions. The first one is 




u = v = - + -J- (3.8) 
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The three other solutions are obtained by the replacement i — > —i, and 
f — > — f. Therefore, it is sufficient to consider explicitly the first case and 
make the necessary replacements at the end for the other solutions. The 
quantity D becomes 

7T 

D = -r== (3.9) 

For this saddle point, / becomes 

/ = i{29 - sin(20)) - 2m (3.10) 

where we have put f = 2cos6>. The fluctuation around this saddle point is 
obtained by the consideration of the second derivatives with respect to u and 
v. They are 

d 2 f d 2 f _ 2i(2 + f)§,2 _ 



du 2 dv 2 \J2 — t t 



d 2 f 4i (2 + f) 



dudv r (2 — f ) i 



(3.11) 



The Gaussian fluctuations around the saddle point produces the inverse of 
the square root of a determinant, which is 

detf" — (—^-) 2 — {—^—) 2 
du 2 dudv 

= 4(2 + f) 4 (3.12) 



Thus we obtain from this result, 

e -iA f (26»-sin26») 



VI - 2«V1 - 2vVDJdetf"N ;i 



S(r) 

~ - - 1 e -iiV(2e- S in2e) (313) 

2 V2isin2^(2 + 2cos#)iV 3 

We now add the other solutions by making the replacements i — > — i and 
f — > — f, which corresponds to 9 — > + 7r. Adding these terms, we have 

o M cos(f - iV(2g - sin2g)) 1 1 
al ] N^V2^29 [ (2 + 2cosfl) + (2-2cos) J l ' J 
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where f = 2cos#; thus 6 = corresponds to r = 2, while 9 = ~ corresponds 
to f = 0. 

The case b) is quite similar. We have 





« = ^ + ^/f^ « = £-^/^ (3-15) 
Using the notation r = 2cos#, we get 

/ = i(26 - sin2#) - iri (3.16) 

and 



The detf becomes 

detf 



d 2 f, 2 , d 2 f 



\2 



dv? dudv 
4(4 -f 2 ) 2 (3.18) 



Adding the four terms obtained by r — > — r and i — > — i, we have 

cos(iV(2g - sin2g) - A% - | ) 
2VriV 3 (4 - r 2 )3 

From these analysis, we obtain the oscillating part of S(t) in the large N 
limit. It is a sum of S a (r) and Sb(r). Noting that the linear part of S(r) in 
r%D is order r/iV 2 = f/N, we find that the oscillating part is a nothing but 



a correction of order ~. However, if we take a derivative, it becomes of the 
same order as the linear term. We also find that when f is close to 2, the 
coefficient of the oscillating part of S(r) becomes large, as shown in ( O J| ) 
and ( |3.19| ), and even diverges at f = 2. Therefore, there should be again a 



crossover near the critical r c = 2N. Up to now we have considered a fixed 
energy E = 0. In the next section, we will take instead an average over E, 
and see that the expression for the form factor simplifies. 



12 



4 Average of the form factor 

We will now consider the average of S(t) over E, by simply integrating over 
E, 



/+oo 
dES(r) (4.1) 
-oo 



Remarkably, we find that this < S(t) > is given analytically in terms of 
known functions. 

From (2.5) and ( |2.18|) , S(t) is written as 



S(t) = J due iU}T p(E, E + lo) 

Thus the integration over E gives simply, 

< S{t) > = JdE J dhe-^+^UoituT) 

= U (-t,t) (4.3) 

Then we write the following contour-integral representation for < 5(r) > 
from ( |7LJ ), 

< S(r) >= / S^(l - £)»(! + £)» ! 



iV 2 J (2vri) 2 v ATit v JVu' (u-w-f) 2 

(4.4) 

Replacing w by r w and u by r v, putting r 2 = x, we have 

< gW >= TP / W («■ - . - *)» (l ~ W (l + ^ (4 ' 5) 

Taking two derivatives with respect to x, we obtain a simple factorized ex- 
pression, 

d 2 < S{x) > e-N j- du e _ ixu _ i N r dv ^ + _z_ N 



cfe 2 N 2 7 2vri AV 7 2vrz v A^ 

'"*[^M^)] a (4-6) 



N 2 L cfe V A^' 
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where Ln(x) is a Laguerre polynomial. Remarkably enough, an identical 
expression has been found earlier, but for a completely different ensemble and 
a different quantity. Indeed the expression fl4.6|) has been found in previous 
work on the Laguerre ensemble of random matrices IJ, in which it was the 
derivative of the density of state. The Laguerre ensemble, also called the 
chiral Gaussian unitary ensemble (chGUE) since eigenvalues appear by pairs 
of opposite signs, has thus a curious relation to the Gaussian unitary ensemble 
(GUE): the form factor < S(r) > in the GUE is related, for any finite N, 
to the density of state p(r) of the chGUE. We have not been able to find a 
direct proof of this exact relation valid for any finite N, without calculating 
both expressions and verifying that they are identical. 

The oscillating behavior of (|4.6|) is similar to that of S(r). The oscillating 
behavior of the density of state for the Laguerre ensemble near the origin can 
be seen in the Fig.2 of (jlTJ. In the large N limit, we know that the oscillations 
of the density of state near zero energy become universal, and are given in 
terms of Bessel functions. In view of the previous correspondence we have 
now to consider the variable r as an energy (although it is a time in the GUE 
problem!). Near zero energy, the density of state of the Laguerre ensemble 
is given by 

p(t) = t[J 2 (t) + J 2 1 (t)] (4.7) 

and it shows an oscillating behavior around one. Consequently, one under- 
stands that the integral of this density of state is proportional to r. This is 
why we have obtained a linear behavior for < S(t) >. However, since the 
density of state for N large is a semi-circle, and not a constant , its integral 
is no longer proportional to r. We have indeed , in the large N limit, by 
integrating the semi-circle line 

<S{r)> = />^l-(^) 2 

T \J 1 ~ (m) 2 arcsin(^) 
" 4iV + 2 (4 ' 8) 

Beyond the critical value r = 2N, it remains equal to one, and it approaches 
smoothly this limit. Therefore, taking into account the fact that the density 
of state is not constant, the singularity is smoothed out. In fig. 3 < S(t) > 
is represented in the large N limit. Near the energy zero, the oscillatory 
behavior of < S(t) >, following from the Bessel functions of (|4.7| ), becomes 
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universal . We have shown in |J, by the same contour-integral representation, 
that there is a crossover from the bulk to the zero energy region, which is 
described universally by the function ((4.7[). We have also found, in a model 
consisting of a lattice of coupled matrices, that this oscillating behavior is 
model-independent. Near r = 2N ', the crossover behavior has been also 
studied in 0. It is given by the square of an Airy function (see eq. (3.37) of 
U); this crossover is also known to be universal. 



5 Time-dependent case 

We now proceed to investigate the time dependent correlation function and 
its Fourier transform, the dynamical form factor. In the large N limit, the 
universal form of this time dependent correlation function has been discussed 



H| [16|, [T7| . We will consider this problem by the contour integral represen- 
tation, which is valid for finite N, and evaluate the form factor S(t) for a 
fixed time t. For a finite t, we will find that S(r) shows different behavior 
compared to the previous linear behavior about r. 

We consider the N by N Hermitian matrix M, which depends upon a time 
t. The time-dependent correlation function is defined by 

p(A, n t) =< ilW(A - M{h))^Tr8{n - M{t 2 )) > (5.1) 

where t = t\ — t 2 and t\ and t 2 are different times. This is written as a Fourier 
transform of the following quantity U(a,/3), 

U(a, (3) = j^< Tre^^Tre^' 2 ) > (5.2) 

We use the new set of variables a and (3 for the Fourier transform variables, 
instead of t\ and t 2 . To avoid the confusion, we use t\ and t 2 as time. 

We show exactly that the correlation function ( |5.1|) reduces to the corre- 
lation function of the two-matrix model in the Gaussian ensemble; the c = 1 
problem is described by the two-matrix model. This correspondence may be 
derived by other arguments [I8|. We here follow the path integral method, 



which can show explicitly that this equivalence to two-matrix model holds 
for any finite N. 
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By considering the following hamiltonian H, 

tf = I T r(p 2 + M 2 ), (5.3) 

where p = M and M is N x N Hermitian matrix, we write U(a, 0) in (|5.2| ) 

as 

U(a,/3) = At < 0\e Htl (Tre iaM )e H{h - h \Tre^ M )e- Hh \0 > . (5.4) 
N 2 

We use the path integral formulation, and we define 



<A\e^ H \B>= [ DMe-^Io {M2+M2)dt (5.5) 

JM(/3)=A,M(0)=B 

Then U(a,(3) is expressed by 

U(a,p) = Jpf dAdB < ^\e Htl \A >< A\{Tre iaM )e H ^- tl) (txe ipM )\B > 

x<B\e- Ht2 \0> (5.6) 

Noting that the ground state energy of the free independent N 2 fermions is 
N 2 /2, we have 

<0\e Hh \A>=e^e-^ A2 (5.7) 
The solution of M = M, becomes 

shf 

M{t) = Bete + ^{A - Bch(3) (5.8) 
sh/3 



Then we are able to write the action in ( |5.5| ) by the matrices A and B, 
2 



-Tr / '(M 2 + M 2 )dt = -Ti(MM)\^ 
! Jo 2 



1 ^[(A 2 + B 2 )chp - 2AB] (5.9) 



2sh/3 

Denoting (3 by a time £, and taking the fluctuation part, we get 



^(«,/3) = ^ 2 / cMdB(Tre^)(Tre^) 

xe -2in TV [( A2 + i?2 ) et - 2 ^ (5.10) 
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Thus the problem reduces exactly to a calculation of the correlation function 
for the two-matrix model, in which matrices A and B are linearly coupled. 

The correlation function for two-matrix model has been studied by D'anna, 
Brezin and Zee ||19|| by the orthogonal polynomial method for finite N. Al- 
though we can use their result, it is more convenient to use the contour 
integral representation for the correlation function. 

By making the change of variables of A, B and a, (3 by a factor ye^ht, 
we obtain a simple expression for ( |5.10| ), 

UfaP) = I f dAdBTre iaA Tre ipB e-^ A2+B2 - 2cAB ^ (5.11) 

Zj J 

where c = e~*. Note that we scaled a and (3 by a factor \J e~*sht, the 
variables A and \x of the two-point correlation function should be modified 
by this factor for the mapping of the time- dependent case to the two-matrix 
model. We now go back to the notation, in which two matrices are given 
by Mi and Mi. We denote the matrices A and B in ( |5.11| ) by M\ and Mi. 
We introduce the external matrix A, which is coupled to matrix M\. The 
Gaussian distribution is given by 

P A (M h M 2 ) = -^e~ H ^ 

#i j2 = ^TrM x 2 + -TrMf - cTrMiM 2 + Ti AM 1 (5.12) 

The density of state p(X) is given by the Fourier transform of 

U A {z) =< ^Tre* 2 ^ > (5.13) 

The calculation of this p(A) is similar to the one matrix case. The integration 
over M 2 , which has eigenvalues is performed by the help of Itzykson-Zuber 
formula. We denote the eigenvalues of Mi by rv The integration over £ 
becomes 

/ dellte - i 3 )e-^~ cN ^ = Ufa - r^Z* (5.14) 

i<j i<j 

Then we are left with the integration about r^, 
1 N r 

U A (z) = — — ^ J drjlfa - r 3 )e-^-c 2 )E^-Nj:^ + i Z r a (515) 



o=l i<j 
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Therefore, by the contour integration, we have by letting a, L goes to zero, 



^-^.(—{\ g ) iv e -7T = ?-^(T^) ( 5 ig) 

We have the same density of state as one matrix case except the scaling 
factor (1 — c 2 ), 

p(A) = v / r^ 2 p (v / r^ 2 A) (5.17) 

where po(X) is the density of state for the one matrix model. In the large N 
limit, this density of state becomes 



iz - -— is*- 



p(A) = y—— y/4-(l-<*)\* (5.18) 



2tt 

which is normalized to be one by the integration. 
The two-level correlation function is given by 

p V)(\ n\ - f f 



(2k) 

where Uo(zi, z<i) is 

C/o(^i, z 2 ) =< I T re l2lMl ^Tre l22M2 > (5.20) 

By the integration over the eigenvalues r, of Mi, and £j of M2, and keeping 
the eigenvalues aj of the external matrix A, we have the following expression, 

T T ( v 1 IL<j( a * - a 3 ~ f - ^>i) - 5v(<*i,a 2 - fy.aa)) 



12^ 

x e 



^ 2 af^2 rii<i( a i-Oi) 

2 2 

l2 l l2 2 c 2 1 2 2 cz l z 1 £ 

T^ aa l^T^ aa 2 ~ 2iV(l-^) ~ 2iV(l-^) ~ jv(l-c^) da l' a 2 (5-21) 



The double sum for cei and a 2 is divided into two parts. The part a\ = a 2 
is written by the contour-integral representation: 



2 

/ i2 2 c \ ' z 2 C2 1 2 2 

^ 1-c 2 ~~ 1-c 2 2AT(l-c 2 ) ~ 2JV(l-c 2 ) ~~ N(l-c' 2 ) (5.22) 



x e 
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The Fourier transform of this quantity becomes by the change of variable u 

to u{z\ + a ), 

I( . i f d Zl dz 2 rdu i N i{zi+Z2c){zi+ ^ )T ^ 

P - N J (2?r)2 f 2m U N J 

z ?+ z 9+ 2cz l z 2 ■ 1 
— = iz\ A— i22M 



x e — 5lvcT^) 1 2 " (5.23) 

This part is further simplified by the change of variables Z\ to 7/j=f (^l — ^2) 
and z 2 to ^7f=f(-22 — czi). By integration over zi, we have 

P{A,N iVfl-c 2 ) 2 J2n?2m [ Nu } 



,2 



x e ^ ^ 2iv(i-c*) (5.24) 
The remaining ai 7^ 02 part is given after letting a 7 go to zero, 

t/oU, = - ^)-(l - " 2 ^ 



ziz 2 / (2vrz) 2V iVw ciW 

2 2 

— "1" _ iz'2 cv z l z 2 

V [1 Ip 1-c 2 1-c 2 2iV(l-c 2 ) 2iV(l-c 2 ) 

[ ciV 2 ( M -,;-Hi)( M -,; + ^) J 

(5.25) 

This expression includes both a disconnected part and a connected part. 
The disconnected part has a factorized form and it corresponds to the first 
term in the bracket. This term, indeed by shifting v — > v/c, becomes the 
product of the density of states p(A) and Therefore, after subtracting 

this disconnected part, we obtain the connected part: 

1 1 dudv izi. Nf iz 2 sN 1 
U (z 1 ,z 2 ) = — r= f 7^-^(1 - -T7-) (1 



iV 2 / (2vri) 2V iW v cAV ( U - V _sa)( w - V + Ma) 

2 2 

iz^u iz^vc 2 i 2 2 

X g 1-c 2 1-c 2 2JV(l-c 2 ) 2JV(l-c 2 ) (5 26) 

where the contour- integrals are taken around u = and v — 0. If we include 
the contour- integration around the pole t> = u — j^, we obtain precisely 
the same term as ( |5.22| ). Therefore we use this representation for the whole 
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expression, including the term of ( |5.21| ), by taking the contour around both 
v = and v = u — iz\/N . 

The expression for the two- matrix connected correlation function p^ (A, p), 
which is obtained by the Fourier transform of Uq(zi, z 2 ), has a factorized form 
when we consider the contribution from u = and v = 0. If z\ and z 2 are 
replaced by Z\ — z\ — iuN and z 2 = z 2 — ivcN, we have a multiplicative form, 



/,),, , / du f dz 2 , Z 2 .at 1 ^id -fiL, — iz 2 a-uNX 

p<f>(X,u) = - <b / — -( — — r 2(1^7 2iv(i-^) 2 ^ 

cN 

dv fdz lfZl , N 1 e -^--JL^ rlZlX ^ c , N ( ^ 27) 



2vriJ 2it K vN' v + % -f 
We write this expression as 

p H (X, p) = -K N (X, p)K N (X, p) (5.28) 
We find the previous first part p / (A, / u) of ( |5.24| ) is also expressed by 

p I (X,p) = ±K N (X,p)e-^ x ~^ 2 (5.29) 

where K^(X,p) is given by 



NXuz 



x e c 5"(i=3) (5.30) 
This kernel Kn(X,p) is written as a sum of Hermite polynomials H n (x), 

K N (X, p) = e-^-^^^ M^m (5.31) 

N ^ c n n! 



where /3 = w|(l — c 2 ). The other kernel K N (X,p) is given 

N nl 

It follows from these expressions that p 7 (A,/i) and p n (X,p) are invariant 
under exchange of A and p. The expressions of ( (5.28Q and ( |5.29| ) for the 
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correlation function by the kernel K^(X,fi) agrees with the result obtained 
by the method of orthogonal polynomials |19 |. The difference is only the 
exponential Gaussian factor in K N (X,fi) and K N (X,fi), and this difference 
disappears for the product of these two kernel. 

In the large N limit, we expect to recover the usual universality. We 
return to the expression of ( |5.26|) and neglect the terms ^ and ^ in the 
denominator. We then exponentiate the powers of N. The integrals over z\ 
and Z2 are Gaussian, and it leads to 

n (Xu) = (1 - c 2 ) r dudv 1 _^ (M+ i^ +(1 - c2)A)2 
P 1 ' W 2nN J (2m) 2 (u-v) 2 

" (wH-i^-Ki-c?)^ 



x e 2 (!-^) 

(1 — c 2 ) r dudv 1 



-N(Mu)+Mv)) 



2ttN J (2ni) 2 (u - v) 



(5.33) 



We use the saddle points of u{X) and u([i), which are the solutions of 
^ = and f 2 - = 0, i. e. 

au av 

w 2 + A(l-c 2 )w + (l-c 2 ) = 

2 1-c 2 1-c 2 . . 

v 2 + fiv + — =— = (5.34) 

c cr 

Taking into account the fluctuations around the saddle point, we have 
ii 1-c 2 1 1 

p (K v) = ~ 2 7t 2 n 2 ^ 77i7~\\ n,( ii\\2 ?m wt ( 5 - 35 ) 



(u(\) - v(n)) 2 fE&EE 

V du 2 dv 2 



where the sum is taken over four different saddle points; note that fi, f% 
vanish at these saddle points. We write down explicitly the expressions for 
u(X) and v(fi) by solving (|5.34|) as 



u(X) = - — - l-X±Jx 2 1 



2 V 1-c 2 



(5.36) 
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where we put A = J and /i — J j^suup. The saddle points become 



l-c 2 r y l-c 

u = iy/T^e ie , -iVT^e- 10 and v = {VT^Fe**, -tyl - c 2 e~^_ Then 
adding these solutions of u and v in terms of 6 and (p, we get from (|5.35| ). 

l-c 2 1 



P n {\n) 



8N 2 7T 2 C COS#COS</3 



2 - (1 + ^)cos(0 - if) 2 + (i + ^) cos (^ + ^) 



v [(l + ^)-fc O s(0-^)] 2 [(1 + ^ + 2008(5 + ^; 

(5.37) 

This expression in the large N limit coincides with the previous result 
jnj. The denominator of ( |5.37| ) does not vanish for A — ► /x. Note that c is 
related to the time t as c = e~*. When t is small, we have c ~ 1 — t. Then, 
the denominator is approximated as (1 — \) 2 + ^ e ~^ ~ t 2 + |(A — /x) 2 , when A 
and /x are small. Note that we have to rescale A and /x for the time-dependent 



case by a factor Ve *sinh£ = y(l — c 2 )/2 as explained in ( |5.11| ). Then we 
have t 2 + |(A — /i) 2 as a denominator and the result agrees with [IS]. In 



the time-dependent case in ( |5.1|) , A and /i are interpreted as one-dimensional 
space coordinates. 

In order to discuss the oscillatory behavior, we return to the expression 
(15.271) . We then change Z2 into Nzc and obtain 

K N (X, „) = cNlp-.f ^-L_ e -*/(-) (5.38) 

J 2m J 2lX U + IZ 

where f(zu) is given by 

c 2 z 2 u 2 

f(z, u) = — + icfxz - Inz + — + Am + lmx (5.39) 

2(1 — <r) 2(1 — +) 

Note that the variables z and w are decoupled, and the saddle point equations 
are simplified. Then, using the previous notaions A = Jj^sxaO and /x = 



Y+^-sin^, we find the relevant saddle points for z and u from the solutions 
of |£ = and If = 0, 

ox ait ' 



+ 1 - c 2 +1 - c 2 ■ 



C C 

u = tVT^e ie , - ivT^AT* (5.40) 
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For the saddle point values, z = ^ 1 c2 e %ip and u = — c 2 e ie , f and the 



c 

,2 



fluctuation determinant o-4|-4 become 



/(z,n) = z(e + ^) + i(e^-e- 2 ^) 
1 A-c 2 \ 



d 2 / 9 2 / y 2c y v/cos^cosy? 

<9« 2 9z 2 

Adding other saddle points values, we obtain 

V /T3^ e f (l-c 2 )(A 2 -^) 

= — ^ V cosftcos^ 
cos[N(h(9) + h(<p)) -(9 + <p)] + ccos[N(h(9) + h(<p))] 

l -±f + cos(9 + <p) 
cos[N(h(9) - h{cp)) -(9- <p)] - ccos[N(h(9) - %>))] 



(5.41) 



1+c 



2c -coa(0-<p) 

(5.42) 

where h{9) = 9 + |sin26' + -^9. Thus we find easily the oscillating behavior 



of Kn(\[Jl) from the integral representation in the large N limit. In |1£ 



this result had been obtained through the integral representation of Hermite 
polynomials. Our derivation is more direct. 

The oscillating behavior of another kernel Kn(\,/a) is obtained similarly 
form (|5.27|) . If we make the smooth average for the product of Kn(X, /i) and 
Kn(X, //), we obtain the previous result (|5.37|) . Also, the expression of ( |5.42j ) 
is the generalization of one matrix result of ( |2.34| ), which may be obtained 
in the limit c — > 1. 

The Fourier transform K(t), the form factor is obtained f rom ( |5.42p . We 
consider the simple case E = 0. 

S(t) = J due iujT p(E, E + u)\ E=0 (5.43) 



Using ( |5.42| ), we have for 9 = 0, simp 



1-c 2 1 



16n*N 2 cosy? [iLjL + sinV] 
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where / is 



/ = [-2cosy?cos[iV h(ip) - y\ + (1 + c 2 )cosNh{(f)} 2 (5.45) 

The Fourier transform of ( |5.43| ) is an integral over u which may be performed 
by taking the residue at the pole to = Then we obtain 

J duoe iuJT K N (0, u)K N (Q, u) ~ re"^ r (5.46) 

Thus we find that the form factor has a linear term in r, but modified by 
e ~ T . Note that the integrand of ( p.43 ) has poles on both sides of the real 



axis the integral does not vanish for r larger than 2a unlike ( |2.7| ). However 
for large r, the form factor becomes exponentially small. The Fourier integral 
of the first term p 1 also becomes exponentially small for large r. Therefore, 
the singularity of the form factor S(t) at the Heisenberg time r = 2N in one 
matrix model is smeared out for the time dependent case. 



6 Correlation functions with an external field 

The contour integral method can be applied to the case in which the Hamil- 
tonian is a sum of a deterministic term and of a random one 0: 

H = H + V (6.1) 

H is a given deterministic term and V is a random matrix with a Gaussian 
distribution P given by 

P(H) = i ( 

= ±_ e -fTr(H 2 ~2H H+Hl) ^ 

Then we are dealing with a Gaussian unitary ensemble modified by a matrix 
source A = —H . 

Therefore, keeping the finite eigenvalues of A and putting H = M, we 
readily obtain the correlation function in the presence of the deterministic 
term A. 
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By the contour-integral method we will show that the connected part of 
the n-point correlation function is also given by the product of the two- point 
kernel Kn(\,[i) in the case of a nonvanishing external source A. 

Let us recall the two-point correlation function, which becomes from 

(PQD 0, 

U A (,ti,t 2 ) = -TT? f /n . A9 e ^iV 2N "- 



N 2 J (2m) 2 ( u - v - f)( u ^ v + 



7=1 



N(u — a^) N(v — a y )' 



By the shift ti — > ti — iwiV, we obtain the kernel of the two point correlation 
function, 



K ^ nc^J ) ^ r -f-*— (6.4, 

The connected part of the two-level correlation function p^'(X, /i) is given by 
— -^Kn(\, fi)K N (fi, A). The Fourier transform of the three-point correlation 
function is given by 

U(h,t 2 , t 3 ) = 4t < Tre itlM Tre it2M Tre it3M > (6.5) 

We evaluate this quantity by the same procedure for the two-point correlation 
function. We integrate out the eigenvalues of M by the use of Itzykson- Zuber 
formula. Then it becomes 



UM ^w^m^tr^ (6 - 6) 



where 6j is given 



°i — a i ^yfi.ai jy°i,a 2 jy0i iCt3 



We have to consider the cases; cui = a 2 = a 3 , oti = a 2 ^ a 3 . These cases 
give the delta- function part as (|2.4j) . 
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When all Oi are different, the contour integral representation is straight- 
forward, 



dudvdw H(u — a 7 — j^)(v — a 7 — ^){w — a 7 v 



(2ni) 3 Y\(u — a 7 )(t> — a y )(w — a 7 ) 



x (u-v)(v-w)(w-u) 1 ^±.±.±. ttl u-u 2V -u 3W (6g) 

(v-w + f) ^lt 2 t 3 



2^-3 i 



Writing the part of the numerators of ( |6.8|) as 



itl ^2n ^1^2 U/ ^2w . ^3, fat 

l(U — V —)(U — V+ —) rrr-J [(V -W T7){ v ~ W + — 



N /y N N 2ny ' N jy ' N J N 2> 

X [( u _ w _ - w + _) _ _] (6.9) 

we have the following disconnected parts, 

U(h)U(t2)U(t 3 ) - U(h)U(t 2 , t 3 ) - U(t 2 )U(tx, t 3 ) - U(t 3 )U(h, t 2 ) (6.10) 

The remaining term is a connected term. We consider the Fourier transform 
of this connected part U c (ti,t2,t 3 ), and make the change of variables, t\ — > 
t\ — iuNfa — > £2 — zfiV and £3 — > £3 — ziuiV. Then the remaining term of 
( |6.9|) becomes simply 

(« + ^ )(v + jfXw + ^) + (u + -Kv + -)(w + — ) (6.11) 

After cancellation of these terms with the corresponding terms in the denom- 
inator, we obtain the connected part of the three point correlation function 

p®(X,p,u), 

pf\X,p, v) = K N (X, p)K N (p, u)K N (u, A) + K N (X, v)K N {v, p)K N (p, A) 

(6.12) 

where K^{X,p) is given by (]6.4f) . Thus we find that the three point corre- 
lation function has the same form as for the unitary Gaussian ensemble in 
terms of the kernel Kn{X,p), but the kernel should be modified as (|6.4j ) for 
the external field case. 
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In the presence of the external source A, the kernel Kn(\,[i) in ( |6.4j ) is 
not expressed as a sum of products of orthogonal polynomials. However, it 
satisfies remarkably the following equation; 



/oo 
d/iK N (X, n)K N {n, v) = K N (\, v) 
-oo 



(6.13) 



The proof of this equation is shown in Appendix B. 

Although the kernel Kn(X, /i) is no longer expressed in terms of Hermite 
polynomials when the external source is present, it still posseses N simple 
eigenfunctions since for n < N — 1, 

K N (X, ^)H n (^N^)e--^ dfx = H n (VNX)e~- x (6.14) 

-oo 

where the H n (x) are the usual Hermite polynomials. However for n > N, 
( |6.14p does not hold. Indeed in the case of zero external source, the right 
hand side of fl6.14|) vanishes. When the external source is non-zero, the result 
is non-zero and is a>i- dependent. The proof of Q6.14 ) is given in Appendix B. 
The n-point correlation R n (Xi, • ■ ■ , A ra ) is defined QTJ] by 



Rn(X, • • • , A n ) 



(N — n)\ 



Pn(Xi, ■ ■ ■ , X N )dX n+1 ■ ■ ■ dX N (6.15) 



/V 



Pn{K;Xn) =<n Tr5 (^- M )> 



(6.16) 



i=i 



Without external source, this n-point correlation function is expressed in 
terms of the kernel K^l^Xi, Xj) as 



where i, j 
has 



R n = det(K N (X l ,X j )) (6.17) 
, n. This result was derived M by the use of (IO) since one 



Pjv(Ai, • • • , A n ) = — det(K N (Xi, Xj)) 



(6.18) 



where i,j = 1,---,N. For a non- vanishing external source, from the rep- 
resentation ( |6.3|) , whose derivation is given in Appendix B, it follows that 
( |6.17|) still holds. This result leads to the universality of the level spacing 
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probability P(s). As shown in a previous work [|7|], the kernel Kn(\, //) has 
a universal short distance behavior, 



sin[iV7r(A - 

K N (\,n)~ * — (6.19) 

[A — {1)71 

The level-spacing probability distribution P(s) is related to the probability 
of having an empty interval with width s E(s), by 

P(s) = ^E(s) (6.20) 

and E(s) is obtained from the integration of Pjv(Ai, • • • , Ajv) in which the 
region — | < Aj < § is vacant. Thus we write 

E ( s ) = IK f°° - I' 2 ) rfA ^(Ai, ■■■An) (6.21) 

Since the kernel Kj^{X,fi) has universal form, same as the case without ex- 
ternal source, we conclude that P(s) becomes same as GUE. 



7 Zeros of the kernel Kn(\, fi) 

Since at short distance the kernel K^(X,fi) exhibits oscillations, and thus 
changes sign, this kernel must have lines of zeros in the (A, /i) plane. Note 
that the solution of the equation 

K N (\,fx) = (7.1) 

are always real in the case of the one matrix model. In the (A, /x) plane the 
solutions of ( |7.1| ) lie on 2N lines as shown in Fig. 2. 

When \i is large enough, -^^(A, \i) is approximated by the product of the 
Hermite polynomial Hn(\) multiplied by // as shown in ( |A.4[ ), and this 
Hermite polynomial has N real zeros. These N real zeros give rise to N non- 
crossing lines in the (A, fi) plane for finite /i. This behavior is a manifestation 
of the short distance universality in the large N limit. The distances between 
these N lines become equal when A — /i is order of 1/N. These N lines are 
parallel to the line A = fi when A and \x are inside the support of the density 
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of state. At the edge, these lines bend, and universality does not hold any 
more. 

We have discussed the time- dependent matrix model, which becomes 
equivalent to the two matrix model. In this two matrix model, the lines of 
zeros of K^(X, n) show a different behavior. In this case, the kernel K^(X, /i) 
is given by ( |5.31| ) or ( |5.32| ). When c ^ 1, the lines of zeros of K N (\,fi) in 
the real (A,/i) plane are not parallel to the A = /i line, as shown in Fig. 3 
on which one can see that some lines turn over. Consequently along the line 
A = — /i, some real solutions are missing. This behavior is related to the fact 
that the short-distance universality does not hold in the two-matrix model. 

In the presence of an external source for the one matrix model, the so- 
lutions of ( |7.1p are not much modified by the source if its eigenvalues are 
smoothly distributed. This is also a manifestation of the short distance uni- 
versality for the oscillatory behavior. When the ths support of the eigenvalues 
of the source are split into separate parts, this is reflected on the position 
of the lines of zeros as shown in Fig. 4. The number of lines of zeros is 
conserved for an arbitrary distribution of the external eigenvalues. 



8 Discussion 

In this paper, we have applied the contour-integral representation for the ker- 
nel which characterizes the correlation functions, for the one matrix model, 
the time-dependent matrix model and in the presence of an external source. 
We have investigated the form factor S(t), which is the Fourier transform of 
the two-level correlation function, by the use of these contour-integral repre- 
sentations. The universality of the two-level correlation p c (A,/i) for A — [i of 
order 1/N, implies immediately the linear behaviour of S(t) in the large N 
limit. We have found explicit deviations from the linear behavior of S(t), and 
found a new surprizing connection to the Laguerre ensemble for the average 
ofS(r). 

Near the Heisenberg time r = r c , a cross-over behavior is observed. For 
the time- dependent matrix model, which we have mapped into an equivalent 
two -matrix model, the universal behavior of the one matrix model is no 
more present and the singularity at r = r c in S(t) is then smeared out. 
This behavior indicates that near the Heisenberg time, the form factor is not 
universal. The non-universality has been pointed out by the authors of M for 
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the case of mesoscopic dirty metals. Our result is consistent with this non- 
universality. Finally we have investigated the zeros of the kernel _K0v(A,/i) 
for the two matrix model, and found differences between the one and two 
matrix models. 

For the matrix model with a non-zero external source, the universality of 
two- level correlation function holds, as shown in a previous paper With 
the technique of contour-integral representations we have also obtained all 
the higher correlation functions. 
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Appendix A: Integral representation for the kernel K^(X,fi) 

It is known that the kernel K]y(X,fi) is a sum of products of orthogonal 
polynomial Pi(X) and Pi(fi) 0). In the case of the Gaussian unitary ensemble, 
the orthogonal polynomials Pi(X) are simply Hermite polynomials. We may 
write these Hermite polynomials as contour integrals 



It is convenient to use at the same time another integral representation for 
these Hermite polynomials obtained by introducing an auxiliary variable t, 




(A.l) 



Their normalization is 




(A.2) 




(A.3) 



30 



The kernel Kn(X, //) is given by 



2n t n 



K N (\,fi) 



(A.4) 



We use two different expressions ( |A.1| ) and ( |A.3| ) of Hermite polynomials in 
the kernel (|A.4|) . The summation of the geometric series give (1— (— ) )/l— ^. 
Then we shift t — > t/y/N and u — > —y/~Nu. 



K n (X,ii) 



00 dt 

-00 27T 



2ni l 



it 



Nu' 



u + — 

N 



e -^-u 2 -^t 2 -it\-Nun 



(A.5) 



This expression coincides with ( |2.22| ) , which has been derived earlier by Kaza- 
kov's method of a vanishing external source. 

This integral representation may also be applied to the two-matrix model. 
The expressions for K N (X, /i) and K N (X, /i) in (|5.27|) are obtained from ( |5.31| ) 
and ( |5.32| ) by the same integral representations of ( |A.1|) and (|A.3| ). 



Appendix B: The properties of the kernel K N (X, /x) in the exter- 
nal source 



We consider the proof of ( |6.13| ). By the integral representation, we write 
the integrand, which is a product of the kernels, by /, 



K N (X,fi)K N (fi, v) 
00 dti / dui f 00 dt 2 
00 2n J 2-ni Joo 2it 
1 1 



dUl JT^ " 7 N N 



<<2 

2m 



Ul 



u 2 



X 



i 2 t 2 

e -^u 2 -^ul-j^-^-it 1 \-it 2 lt-Nu 1 ^-Nu 2 v 



(B.l) 



Making the shift £2 — ► ^2 + iNu%, and integrating J over //, we obtain the 
5(t 2 ) function. Thus the integral of I becomes 



dfil 



-00 2tt / 2vri ' °— ' ^ 



x e 2 ^ 



^ — iti A— u| — Nu2 v 



2ni \ Ly u 2 - a 7 (u x + ^)(u 2 - «i) 

(B.2) 
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The contour-integral of U\ can be performed around U\ = — ^ since U\ 
appears only in the denominator. Then, we obtain 



= K N (\u) (B.M) 



' u 2 - a 7 u 2 + ^ 



The integral equation of ( |6.14| ) is verified similarly. We evaluate first the 
following integral involving Hermite polynomials, 



/OO f* OO 

H n (VNfi)e-^ +u)2 dfi = / H n {fi)e-^ +y ^ u)2 
-oo J—oo 



du 



(-VNu)\ _i u 2 dfi 



i n — 1 

= V2^(-u) n N— (B.4) 

Using the expression of ( |6.4| ) for K^(X,fi) with ( B.4 ), we obtain 

r i — n 9 
/ K N (\,fi)H n (VNfi)e~-^ dfj, 



v ; i-oo 2tt J 2ni ^ u-a 1 J u+§ v ; 

When n < N, this contour-integration of u converges for \u\ — > oo; we 
may then take the residues of the poles outside of the contour , i.e. the pole 
u = — £ instead of the a 7 's. Then, by evauating the residue, and performing 
the t integration, we obtain 

J—oo J—oo v Z7T iV 



00 dt f du + % A-u 



= H n (VN\)e~% x2 (B.6) 

Note that the contour-integration around infinity does not converges for n > 
N; we cannot take the poles outside of the contour any more and we obtain 
a different result. 

The kernel Kn(\,/i) for a non-zero external source can be written as 
a determinant; this is useful for numerical calculations. Since the kernel 
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Kn(\,h) satisfies ( |B.6| ), we write an expression for Kn(X,/j) as a determi- 
nant, in which the variable \x appears only in the first row of the N x N 
matrix. The eigenvalues of external source cij appears only in the i-th col- 
umn of the matrix. This is related to the fact that the exchange between aj 
and a,j in ( P-7|) does not affect K^{X,p). We can write the matrix element 
of the i-th row as a polynomial of order of (i - 1) for A. The first three rows, 
for example, are expressed by 



K N (\,ji) 



N i2 
e -T A 



27T Ui<j{ai - ajj 

I e -^a\-Na lt i 



x det 



\ 



A + a\ A + a 2 

(A + ai ) 2 + ^ (A + a 2 ) 2 + ^ 



^a? N -Na N Li \ 



A + Oat 

(A + a N ) 2 + i 



(B.7) 



The matrix element is (A + a,,)* -1 when i is even. When i is odd, then 
rriij is given by (A + a^)* -1 + C, where C is a constant depending upon N. 
This constant is determined to be consistent with ( JB.6| ). 
As a simple example for N=2, we obtain 



'2-k (oi — a 2 ) 



~ x2 det 



A + ai 



A + a 2 



(B.8) 



It may be interesting to note that there exists values of the cij for which 
the two level correlation function p c (X, //) = — K N (\, fi)K N (fi, A) , which is 
normally negative, may become positive. For instance, in the case N = 2, 



we consider \i - 
-l/v / 2T(Ae- a ? 



0, a 2 = 0. In this case, we obtain from ( [B.8|) , K 2 (\,0) 
- A — ai)e~ A ja\ and ^(0, A) = \ j\p2/K. Then, we have 



Pc(A,0) 



87ra! 



(Ae ai - A - ai)e 



(B.9) 



This expression may take positive values, for example when a% — 1 and 
A < —1.6. The expression (p.7| ) is useful for locating numerically the zeros 
of K N (X,fi). 

Finally let us mention that the contour-integral technique which we have 
used here may be extended to the n-point functions. If we define 
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jj,. . s _ f II duj SL (uj a 7 ) -j-j- Uj n + n) 

x J](^-^)-^ e -^S*?-iE*^ (B.IO) 



then -R n (Ai, • • • .A n ) is obtained by Fourier transform of U(ti, • • • , t„) as 

i?n(Ai,---,A„) = Ju(t 1 ,---,t n )e-^ t ^dt l (B.ll) 
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Figure caption 

• Fig. la, A quasi-linear behavior of S(t), N = 16 

• Fig. lb, The derivative of S(r) of Fig. la. 

• Fig. 2. The zeros of K N (\,/i) are plotted as lines in the real (A, //)- 
plane, for N = 5. 

• Fig. 3. The lines of zeros of K N (\, //), N = 5 for the two-matrix model 
with c—\. 

• Fig. 4. The lines of zeros of K N (\,fi), N = 5 for the external source, 
di = —2, a 2 = 2, a 3 = 2.25, a 4 = 2.5 and a 5 = 2.75 
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